Spatio-temporal vortex beams and angular momentum 
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We present a space-time generalization of the known spatial (monochromatic) wave vortex beams 
carrying intrinsic orbital angular momentum (OAM) along the propagation direction. Generic 
spatio-temporal vortex beams are polychromatic and can carry intrinsic OAM at an arbitrary an- 
gle to the mean momentum. Applying either (i) a transverse wave-vector shift or (ii) a Lorentz 
boost to a monochromatic Bessel beam, we construct a family of either (i) time-diffracting or (ii) 
non-diffracting spatio-temporal Bessel beams, which are exact solutions of the Klein-Gordon wave 
equations. The proposed spatio-temporal OAM states are able to describe either photon or electron 
vortex states (both relativistic and nonrelativistic) , and can find applications in particle collisions, 
optics of moving media, quantum communications, and astrophysics. 

PACS numbers: 42.50.Tx, 03.65.Pm 



I. INTRODUCTION 

Wavefront dislocations (phase singularities or optical 
vortices) were introduced in 1974 in a seminal paper by 
Nye and Berry [1] which gave rise to the field of singu- 
lar optics [2]. A non-singular wave locally looks like a 
segment of a plane wave with parallel wave fronts, well- 
defined phase, and current flowing in the direction of 
wave propagation. Nye and Berry demonstrated that, in 
general, the wavefronts in 3D space contain topological 
dislocation lines, akin to those in crystal lattices, where 
the phase becomes singular and currents coil around the 
lines forming optical vortices. The dislocations are char- 
acterized by a quantized vortex strength, also called topo- 
logical charge. One can distinguish screw and edge wave- 
front dislocations where the vortex lines are, respectively, 
parallel and orthogonal to the wave propagation direc- 
tion. In general, mixed edge-screw dislocations, with the 
vortex core tilted with respect to the phase fronts, are 
possible. Importantly, only screw dislocations are generic 
in monochromatic wave fields, while generic edge-screw 
dislocations require additional temporal variations [1,2]. 

In 1992 Allen et at revealed a close relation be- 
tween optical vortices and wave orbital angular momen- 
tum (OAM) [3]. They showed that axially-symmetric 
higher-order optical beams in free space bear screw dis- 
locations along their axes and possess intrinsic OAM di- 
rected along the beam axis. This OAM is quantized with 
its magnitude (in units of h per one photon) being equal 
to the vortex charge. Thus, monochromatic vortex beams 
represent the OAM states of light, and nowadays they 
play an important role in classical and quantum optics 
[4]. Remarkably, the above phase-dislocation-, vortex- 
, and OAM- properties are generic for all types of lin- 
ear waves, both classical and quantum, independently of 
their nature. In particular, vortex beams with OAM were 
recently described and generated in electron microscopes 
[5,6], as well as employed in acoustics [7]. 

For monochromatic vortex wave beams, the intrinsic 
OAM is collinear to the momentum, with its projection 



on the beam axis (helicity) being quantized. In other 
words, such OAM behaves very similar to the spin of a 
massless particle (even for nonrelativistic electrons [5]). 
Then, a natural question arises: 11 Can a generic wave 
packet or beam carry a well-defined intrinsic OAM in an 
arbitrary direction, i.e., tilted with respect to the propaga- 
tion direction?". It would seem that the direction of the 
intrinsic OAM is associated with the direction of the vor- 
tex line, and the tilted OAM should appear in states with 
tilted vortices, i.e., mixed edge-screw wavefront disloca- 
tions. As we pointed out, such dislocations are generic for 
polychromatic fields and our question essentially requires 
studying vortex beams and OAM in space-time [8,9]. 

In this paper, we address the above question and ex- 
tend the concepts of intrinsic OAM and vortex beams 
to polychromatic states in space-time. We analyze the 
scalar Klein-Gordon wave equation, so that its relativis- 
tic and nonrelativistic limits describe both massless op- 
tical fields and massive Schrddinger particles, assum- 
ing that the polarization effects can be neglected. We 
show that the existence of spatio-temporal vortex beams 
with tilted OAM follows from the requirement of rel- 
ativistic invariance and can be obtained via Lorentz 
transformations of the usual spatial vortex beams [9]. 
Hence, the spatio-temporal OAM wave states can nat- 
urally appear from moving sources emitting stationary 
(monochromatic) vortex states in their rest frames. As 
such, our results have implications in optics of moving 
media [10], quantum communications with satellites [11], 
collisions of high-energy particles with OAM [12], and as- 
trophysical applications of OAM of light [13,14]. 



II. MONOCHROMATIC BESSEL BEAMS 

We start with the Klein-Gordon wave equation in units 
h = c=l: 



-df + V 2 - m 2 ) ip = 



(1) 



2 




where i/}(t, r) is the scalar wave function. The plane- wave 
solutions of this equation are: 



tp oc exp [i (— uj t + k • r)] = exp (ifc^r M ) , 
cu 2 - k 2 = k^K = m 2 , 



(2) 



with fc M = (cJ,k) and r M = (t, r) being the standard 
four-vectors in the Minkowski space-time with signature 
(—,+,+,+). The second equality in Eq. (2) is the dis- 
persion relation which determines the mass hyperboloid 
(or the light cone at m = 0) - a hypersurface in momen- 
tum & M -space, where we only consider the positive-energy 
domain ui > 0. We deal with the Klein- Gordon equation 
because it possesses relativistic space-time symmetry and 
is able to describe both quantum massive particles of dif- 
ferent energies and classical waves. (Indeed, for k -C to 
it can be reduced to the Schrodinger equation, whereas 
in the relativistic limit k 3> to it becomes the usual wave 
equation). Alongside the wave function i/j(t,r), we will 
use its plane-wave (Fourier) spectrum ip(u>,)t), defined 
with the explicit delta- function of the dispersion (2): 



4> (r») <x H> (k^) <5(fc% - to 2 ) e ik " r >'d 4 k>* 



(3) 



The simplest monochromatic solutions of wave equa- 
tions, which carry OAM, are Bessel beams [15]. These 
represent eigenmodes of the z-components of the momen- 
tum, p = — iV, and OAM, L = r x p, and can be con- 
structed as a superposition of multiple plane waves (2) 
with fixed uj = ujq and k z — k Z Q. The fixed frequency im- 
plies that the wave vectors form an isofrequency sphere 
k = ,/uj 2 



to = fco in k-space, whereas a fixed k z cuts 
fcj_o on this sphere, see 



k 2 



a circle k± — y^ui 2 — m 
Fig. la. Hereafter (k±,cj),k z ) denote the cylindrical co- 
ordinates in k-space. The relative phases of plane waves 
can grow around the circle and form a vortex exp (i£(f>) 
with phase increment 2ir£ around the loop, whereas the 
vortex charge £ = 0,±1,±2,... determines the order of 
the Bessel beam and its OAM. Thus, the plane-wave 
spectrum of the i'th-order beam can be written as 

ipe (uj, k) oc S(k — ko) 5(k± — k± ) exp(i£0) . (4) 

The corresponding real-space field (3) yields 

tpz (t,r) cx J w (k± r±)exp[i(£(p + k z0 z-uj t)] , (5) 

where J n (C) is the Bessel function of the first kind, and 
(rj_, tp, z) are cylindrical coordinates in real space. In the 
vicinity of the beam axis, k±or± <C 1, 

ipe oc [x + isgn(f) y] W exp [i (k z0 z - w t)] , 

which demonstrates a screw wavefront dislocation of 
strength t on the axis [1,2]. 

Bessel beams (5) are non-diffracting, i.e., \ipe\ does 



not vary with z and t [15] 
modes of the operators p z 



Evidently, these are eigen- 
— id z and L z — —id v , with 
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FIG. 1: (color online), (a) Plane-wave spectrum (4) of a 
monochromatic Bessel beam with t — 2: a circle on the isofre- 
quency sphere k — \J uj 2 , — m 2 in k-space. Relative phases of 
the constituent waves form a charge-2 optical vortex encoded 
by colors, (b) Transverse real-space distributions of the in- 



tensity p 



current j = Im(V'Vi/>) = pV$, and phase 



3> = argi/i for the Bessel beam (5) corresponding to Fig. la. 
The phase singularity on the beam axis represents a screw 
wavefront dislocation accompanied by a vortex of the cur- 
rent. The parameters used here are uio/m = 2, k±o/m — 1, 
k z o/m — v2, and the dimensionless coordinates are x = k±ox 
and y = k± y. The normalized momentum and OAM, Eqs. 
(6) and (7), integrated over the visible area, yield p = k z oe z 
and L = £e, . 



corresponding eigenvalues k z and £. Figure lb shows the 
transverse spatial distributions of the intensity (proba- 
bility density) p = \ip\ , current (momentum density) 
j = Im(^*V-0) = pV$, and phase $ = arg^ for the 
beam (5). It is the screw phase dislocation at r± = and 
the accompanying vortex current that produce a non-zero 
OAM of the field. 

To characterize the angular momentum, it is impor- 
tant to also describe the expectation values of energy- 
momentum and coordinates of the beam [9,16]. The nor- 
malized (per one particle) expectation values of the en- 
ergy, transverse coordinates, momentum, and OAM for 
the Bessel beam are: 



E = 



(-01 -0) 
(-01 -0> 



U>Q. 



k z oe z , L 



= 0, 



£e z .(6) 



Hereafter e a denotes the unit basis vector of the corre- 
sponding a-axis, whereas the subscript "_L" indicates the 
transverse (x, y)-components of a vector. As expected, 
the OAM is collinear with the momentum, L || p. The in- 
ner product in Eq. (6) implies integration over the proper 
area of space. Formally, the Bessel beams are not local- 
ized, i.e., cannot be normalized in the whole space. How- 



3 




nil 



phase: 

71 



3rc/2 



(b)io 



-10 





p j 

o 

t=0 


t=15 









-10 



10-10 



~ 5 



10-10 



-5 0-5 10 



FIG. 2: (color online), (a) Plane- wave spectrum of a polychromatic Bessel beam with I = 2 and spectrum (9) shifted by 
Ak e x in k-space. Here the semitransparent circle represents the initial monochromatic beam, and the isofrequency sphere 
has radius y&JJ + Afc 2 . (b) Transverse real-space distributions of the intensity p, current j, and phase $ for the beam (a) at 
different times t. The beam moves in the positive ^-direction and experiences time-diffraction deformations. The fork-like phase 
patterns in the beam center represent moving edge-screw wavefront dislocations. The parameters are k z o/m — 1, k±o/ m = 1-4, 
Ak/m = 2, u c /m = 2.45, and we use the dimensionless variables x = k±ox, y = k±oy, and t = 0J c t. The normalized expectation 
values, Eqs. (6) and (7), numerically integrated over the same area around the moving beam center yield Eqs. (11) with a good 
accuracy: Ri ~ (Ak/uj c ) te x , p ~ k z oe z + Afce^, L ~ I e z ]\ p. 



ever, this can be overcome by substituting delta-functions 
in the spectrum (4) with arbitrarily narrow Gaussian ex- 
ponents, and the normalized expectation values (6) con- 
verge to finite values. An alternative form of Eqs. (6) 
using the probability density p and momentum density j 
distributions can be written as [17] 



R i - J . = 0, p = Y „ - = k z0 e z , 



Jpd 2 r ± 
L= /(rxj)dV =fez _ 
J pa z r± 



J pd 2 r± 



(7) 



This makes it clear that it is the circulation of the current 
j (shown in Fig. lb) that produces the intrinsic OAM 
of the beam. Note that the proper probability density 
and momentum density for the relativistic Klein-Gordon 
equation differ from the "naive" p and j used here [9,18]. 
However, this does not affect our consideration and will 
be discussed below in Section V. 

In practice, the Bessel-beam spectrum (4) can be sim- 
plified and well approximated by a finite number N 3> 1 
of plane waves distributed over the spectral circle in k- 
space (see Fig. la) [14]. In this manner, 



ipi (t, r) cx 



1 N 

—j= exp [i {—u) n t + k„ • r - 

v ™ n=l 



■e<M], 



u n = wo, k„ = (fcj_ cos <f) n ,k± sin <f> n , k z0 ) , 



(8) 



The discretization (8) is important for interferometric ap- 
plications [14] and approximates well the beams (5) in 
the restricted area around the axis: e.g., k±Qr± < 10 for 
N = 20. In Figure 1 and throughout this paper we use 
the superposition (8) of N = 30 plane waves for numer- 
ical simulations of Bessel beams, because in some of the 
cases described below the beams cannot be characterized 
analytically. We verified numerically that the normalized 
energy, momentum, and OAM of such beam are in perfect 
agreement with Eqs. (6) and (7) when the integration is 
performed over the area S = {\k±ox\ < 10, |fc_i_oj/| < 10}- 



III. SPATIO-TEMPORAL BESSEL BEAMS 
WITH SHIFTED SPECTRUM 

From the above picture of the monochromatic Bessel 
beams, one can see that the expectation values of the 
OAM and momentum are collinear not by chance. In- 
deed, assuming a circular plane- wave distribution in mo- 
mentum space (Fig. la), geometrically, p represents the 
radius- vector of the center of the circle, whereas L points 
in the direction normal to the circle and has a magnitude 
equal to the vortex charge I. Obviously, the collinear- 
ity L || p holds true for any circle on the isofrequency 
sphere, i.e., for any axially-symmetric monochromatic 
beam. However, as soon as we abandon the monochro- 
maticity constraint, it is possible to construct a Bessel- 
beam-type solution with L \ p. Indeed, any circle in k- 
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FIG. 3: (color online). Spatio-temporal representation of the transversely-moving Bessel beam from Figure 2. (a) Here its 
spectrum is plotted on the mass hyperboloid in the (u),k x , fe !/ )-space where it is seen that it also forms a vortex loop if projected 
onto the (cj,k v ) plane, (b) Real-space evolution shown here in the form of the distributions in the (t,y) plane, for different 
values of x, cf. Fig. 2b. The corresponding components of the four-current j' M are (j°,j y ), where j° = —Im(t/j*dtip)- 



space can serve as a plane-wave spectrum for the beam. 
The only fundamental constraint of the mass hyperboloid 
(2) can always be satisfied by choosing the corresponding 
frequencies uj (k) = y/m 2 + fe 2 . Thus, considering differ- 
ent spectral circles with vortex phases (see an example in 
Fig. 2a), we obtain different vortex solutions of the wave 
equation (1). It is easy to show that the center of the cir- 
cle and the normal to the circle still represent the mean 
momentum and intrinsic OAM of the beam, respectively. 

The simplest transformation shifting the spectral circle 
away from the isofrequency sphere is a uniform transverse 
shift (e.g., along the x-axis) of all the wave vectors in the 
Bessel-beam spectrum (Fig. 2a): 



k^k + Ake x , w -> w(k) = y + Afc 2 + 2k x Ak. 

(9) 

Solutions with spectrum (9) are polychromatic because 
of the u{k x ) dependence (here k x is the wave vector com- 
ponent of the original monochromatic beam). Remark- 
ably, the beam spectrum (9) represents a circle with 
an azimuthal phase gradient not only when projected 
onto the (k x ,k y ) plane, but also a similar vortex loop 
around u> — uj c = \JuJq + Ak 2 in the (w, k y ) plane. Thus, 
this is a vortex loop on the mass hyperboloid in (w,k) 
space, which demonstrates its spatio-temporal nature, see 
Fig. 3a. 

Formally, the wave function of the beam can be written 
as a Fourier integral (3) with the delta-function spec- 
trum (4) shifted by the transformation (9). However, 
such integral cannot be evaluated analytically in the gen- 
eral case. To understand properties of the beams with 
a shifted spectrum, we performed numerical simulations 
using the discretization (8) with parameters correspond- 



ing to Eq. (9): 

U! n = a/cJq + Afc 2 + 2fcj_o Ak COS (j) n , 

k n = (Afc + k±Q cos <fi n , k±Q sin <j> n , fc z o) ■ (10) 

The results are presented in Figure 2b. We see a 
Bessel-beam solution which is still homogeneous in the 
z-direction [k z — k z o) but evolving in time. At t = 0, 
the intensity distribution represents the Bessel function: 
p(r,t = 0) oc \Ji£ | (fc±o r ±)] ! but the current and phase 
distributions differ significantly as compared with the 
monochromatic beam (5) (Fig. lb). Indeed, there is 
a net transverse current flowing in the x-direction (be- 
cause of the mean (k x ) = Ak), whereas the phase dislo- 
cation in the beam centre becomes of mixed edge-screw 
type. This is seen from the fork-like pattern of the phase 
fronts instead of the radial pattern characteristic of screw 
dislocations. Calculating the wave-field distributions at 
( ^ 0, we see that the beam moves in the transverse 
x- direction and experiences shape distortions. The de- 
formations represent diffraction in time caused by the 
accumulated phase difference between waves with differ- 
ent frequencies (10). In particular, Fig. 2b shows that 
this diffraction breaks the £th-order phase dislocation at 
t = into \£\ basic dislocations of strengths at t ^ 0. 
This can be regarded as a temporal manifestation of the 
general instability of higher-order vortices with respect 
to perturbations [1,2,19]. 

The moving edge-screw dislocation reveals a spatio- 
temporal character of the shifted-spectrum vortex beam. 
To illustrate this, we plot the beam density, current, and 
phase distributions not only in the (x, y) plane for differ- 
ent values of t (Fig. 2b), but also in the (t,y) plane for 
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FIG. 4: (color online), (a) Plane- wave spectrum (15) of a Lorentz-boosted Bessel beam with I — 2 in k-space. Here the 
semitransparent circle represents the initial monochromatic beam, whereas the isofrequency sphere has radius y k^ + (■yvujo) 2 . 
(b) Transverse real-space distributions of the probability density p, current j, and phase $ for the beam (14) corresponding 
to (a), at different times t. The beam moves uniformly in the positive x-direction without diffraction. The fork-like phase 
pattern represents a moving edge-screw wavefront dislocation. The parameters are uio/m = 2, Ko/m = 1, Ak/m = 2, v = 0.8 
(7 = 5/3), and we use the dimensionless variables x = k±o x , V = k±oy, and t — cj c t. The numerically-calculated normalized 
expectation values, Eqs. (6) and (7), correspond to Eq. (18) with good accuracy. 



different values of x, in Fig. 3b. One can see quite sim- 
ilar "moving" beams with edge-screw phase dislocation 
and circulating probability current in both the (x, y) and 
(t, y) planes. Thus, in general, this is a vortex "hyper- 
beam" in Minkowski (t, r) space-time. 

The transverse motion of the beam implies a non-zero 
transverse momentum, p x ^ 0, whereas the OAM is still 
determined by the z-direction of the vortex axis (which 
also moves along x with time). Calculating the expec- 
tation values (6) in the momentum (Fourier) represen- 
ta tion and a ssuming paraxiality of the beam, fc^o ^ 
yj k 2 + Ak 2 , we obtain 



IV. LORENTZ-BOOSTED SPATIO-TEMPORAL 
BESSEL BEAMS 

It turns out that it is possible to construct a family 
of spatio-temporal beams which are free of the tempo- 
ral diffraction and allow a simple analytic description. 
Note that Lorentz transformation to a moving reference 



frame, 



A M „(v) r" (v is velocity), provides a boost of 



the four- momentum in the spectrum, fc M — > A^ u (— v) k v ' , 
and can generate polychromatic moving solutions. Since 
the Klein-Gordon wave equation is Lorentz- invariant, the 
Lorentz-boosted Bessel beams (4) and (5) will also rep- 
resent exact solutions of Eq. (1): 



E 
P 



)l + Ak 2 , Rj 



Afc 



k zQ e z + Ak e x 



(11) 



Here the beam centroid moves in the transverse di- 
rection according to the free-space equation of motion: 
dR±/dt = p±/E. Numerical calculations of the expec- 
tation values in the representation (7), performed using 
the same area around the instantaneous transverse beam 
centroid Rj_, showed good agreement with Eq. (11) even 
for the non-paraxial parameters used in Fig. 2. Thus, 
Figures 2 and 3, together with Eqs. (11), demonstrate 
the existence of spatio-temporal vortex beams with non- 
collinear momentum and intrinsic OAM: L l/f p. 



$ (r") = ^ [A^(v) r»\ , # (k v ) = fa |A"„(-v) k" 

(12) 

Equations (12), with Eqs. (4) and (5), offer a new family 
of moving vortex beams parameterized by the velocity v. 

It is easy to see that the Lorentz boosts in the longi- 
tudinal z-direction keep the monochromatic Bessel-beam 
family (4) and (5) invariant and only transform their pa- 
rameters wo, k Z Q, and k±Q. Therefore, we consider only 
the non-trivial case of the transverse boost, say, in the 
x-direction. Assuming v =ve x , we have 



A"„(v) 



7 — 1>7 
-wy 7 
10 
1 



(13) 
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FIG. 5: (color online). Spatio-temporal representation of the transversely-moving Bessel beam from Figure 4. (a) Here its 
spectrum is plotted on the mass hyperboloid in the (w, k x , fcjj-space where it is seen that it also forms an elliptical vortex 
loop if projected onto the (ui,k y ) plane, (b) The real-space evolution of the beam is shown here in the form of the density, 
current, and phase distributions in the (t, y) plane for different values of x, cf. Fig. 4b. The corresponding components of the 
four-current j M are (J°,jy), where j° = —Im(ip*dtip)- 



with 7 = 1/vl — v 2 being the Lorentz factor. Explicitly, 
the moving Bessel beam (12) with Eq. (13) is obtained 
via the substitution t — > j(t — vx), x — > 7 (x — vt) in 
Eq. (5), and it takes the form 

ip' e (t,r) oc J\ e \(kxor'x)exp[i(£<p' + k zQ z - uj t')} , 
r± = \A 2 ~ vtf + y 2 , <p' = tan -1 > 

t'=j(t-vx). (14) 

At the same time, the Lorentz transformation (12) and 
(13) of the spectral characteristics yields [cf. Eq. (9)]: 

k x ^ (k x + vujq) , loq -> uj (k) = 7 (uj + vk x ) . (15) 

In k-space, this represents a shift of the spectral circle 
by the distance Afc^ = 7VW0 and a stretch with factor 
7 along the fc^-dimension, see Fig. 4a. This stretch- 
ing of the spectrum ensures the linearity of the u(k x )- 
dependence which causes no diffraction in time. Plotting 
the Lorentz-boosted spectrum (15) in the (u>, k) space, 
one can see an elliptic planar loop on the mass hyper- 
boloid, Fig. 5a. One can consider this as a vortex loop 
both when projected onto the plane and also on the 
(k x ,k y ) plane, which demonstrates the spatio-temporal 
character of this vortex. The spectral discretization (8) 
for the Lorentz-boosted beam acquires the following form 
according to Eq. (15): 

uj n = 7(^0 + vkj_ cos (j) n ) , 

k„ = (7(fcio cos <f> n + vuq) . k ±0 sin tj> n , k z0 ) . (16) 

Figure 4b shows the real-space distributions of the in- 
tensity, current, and phase in the Lorentz-boosted Bessel 



beam (14). This beam moves in the transverse in- 
direction with velocity v = ve x and has an elliptical 
profile (stretching of k x yields the Lorentz contraction of 
the s-dimension with factor 7). The wave function (14) 
near the beam axis, fcj_o r l <C 1, is approximated by 

tp'i oc [7(2; — vt) + i sgn(£) y ]'^' 

x exp[i (k z0 z + o; 7 (vx - t))] . (17) 

It has the form of a mixed edge-screw dislocation [1,2], 
i.e., a spatio-temporal vortex in the (t,x,y) space, with 
its singularity line {y = 0, x = vt} being parallel to the 
z-axis and moving along the a;-axis. Figure 5b illustrates 
the spatio-temporal character of the Lorentz-boosted 
beam (14) by showing its evolution in the (t,y) plane 
for different values of x. One can see that the x-moving 
vortex appears in the (t, y) plane entirely similar to the 
t-moving vortex in the (x, y) plane. 

Similar to the example discussed in the previous Sec- 
tion, the transverse current and motion of the Lorentz- 
boosted beam implies a tilted momentum of the field, 
p, x 7^ 0, whereas the OAM is still determined by the in- 
direction of the vortex axis. Calculating the expectation 
values (6) and (7) and assuming paraxiality of the beam, 

-fk± <C y k 2 M + (jvujq) 2 , we obtain: 



E ~ lo c = \Ju 2 + (7WW0) =7 w , R_L~uie x , 

p ~ k zQ e z + 7wa;oe x , L ~ ie z . (18) 

Here the values of the transverse coordinates and mo- 
mentum correspond to the linear motion of the beam, 
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FIG. 6: (color online). Same evolution of the Lorentz-boosted Bessel beam (14) as in Fig. 4, but here with the scaling ip — > y/uiip 
of the spectral Fourier amplitudes. Thus, p = \ij)\ 2 becomes the true relativistic probability density, and j = Im (i/>*Vi/>) becomes 
the energy current (momentum density). This causes the relativistic Hall effect [9]: a tiny y-deformation and subwavelength 
shift (19) of the beam which is seen here owing to the nonparaxial character (small radius) of the beam. The dot indicates the 
centre of gravity of the beam according to Eqs. (18) and (19): Rj_ — (vt, —v(./2uio), and the small ^-deviations of the dot are 
due to the non-paraxial corrections. 



dH±/dt = p±/E, whereas the transformation of the 
OAM, which acquires the factor (7 + 7 -1 ) /2, owes its 
origin to the elliptical relativistic deformation of the 
beam [9,20]. Indeed, the OAM is given by the cross- 
product of the position and momentum, L z — xp y — yp x , 
and both summands xp y and — yp x provide equal contri- 
butions of L z /2 to the intrinsic OAM of the cylindrical 
beam (5). In the transversely-moving beam (14) and 
(15), the Lorentz transformation results in the contrac- 
tion of the coordinate x and expansion of the momentum 
component k x (Fig. 3), and the above two summands ac- 
quire factors 7 _1 and 7, respectively. Thus, Eqs. (12)- 
(18) and Figures 4 and 5 demonstrate the existence of 
spatio-temporal non- diffracting Bessel beams with non- 
collinear momentum and intrinsic OAM, L !/( p. 



V. RELATIVISTIC HALL EFFECT 

It should be noticed that the last equation (18) con- 
tradicts the Lorentz transformation of the angular mo- 
mentum L of a point particle. Indeed, the Lorentz trans- 
formation of the original OAM L = £e z should yield 
L' = j£e z after the boost (13). This apparent paradox 
is considered in detail in [9]. It is resolved when one 
takes into account the "relativistic Hall effect", i.e., the 
transverse y-deformation and shift of the OAM-carrying 
object observed in an a;-moving frame. It turns out that 
the geometric center of such object undergoes a trans- 
verse shift of its centroid: 



we obtain 



Y 



v vi 
2Eq 2cjq 



(19) 



Together with the momentum component p x ~ ^vujq, 
Eq. (18), the shift (19) generates extrinsic OAM L (cxt ) = 
-Yp x = jv 2 e/2. Adding it to the intrinsic OAM (18), 



L (int) + L (oxt) _ 7 + 7 1 1 



7i> 



(20) 



in agreement with the Lorentz transformation. 

We did not observe the shift (19) in Figures 4 and 
5, because we used the "naive" density p — \if)\ and 
current j = Im (if*'Vifj). Indeed, p does not describe the 
density of particles in the case of the relativistic Klein- 
Gordon equation (1) because the volume integral J pdV 
is not Lorentz- invariant. The proper particle density and 
current are given by [18] 



Pp 



-Im(^*atV0, jp=Im(^*VV0: 



(21) 



whereas the energy density and current (momentum den- 
sity) are: 



1 r 



PE 



|cW| 2 + |W| 2 +/i 2 M 2 



j E = -Re[(a 4 v>)*(W)] 



(22) 



For plane waves, this yields simple w-scalings: pe = p, 
]e = wj, and pp = to p, which make no difference for 
monochromatic beams in the rest frame. However, a 
Lorentz transformation to the moving frame affects these 
distributions via local variations of the frequency. In 
other words, one can use the "naive" density p and cur- 
rent j, instead of pp and jg which should be used in 
Eqs. (6), but the amplitudes of the Fourier spectral com- 
ponents should be scaled as ip — > y/to if. 

In Figure 6 we plot the spatio-temporal evolution of 
the Lorentz-boosted beam with the Fourier amplitude 
scaling ip — > *JZi tp taking into account the relativistic 
Hall effect [9]. The y-deformation and shift (19) are 
clearly seen owing to the chosen nonparaxial relativistic 
parameters of the beam. However, for paraxial or nonrel- 
ativistic beams the tiny shift (19) (less than a fraction of 
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the wavelength) is practically imperceptible. Moreover, 
we emphasize that both the spatio-temporal beams with 
scaling tp — > y/ujtp and those without it are exact solu- 
tions of the Klein-Gordon wave equation and both can 
exist. 



VI. CONCLUDING REMARKS 

We have demonstrated the existence of spatio- 
temporal vortex wave beams in free space. Such beams 
move uniformly in the transverse direction and carry in- 
trinsic orbital angular momentum non-collinear to their 
mean momentum. Furthermore, the spatio-temporal vor- 
tex around the beam axis appears as a moving edge-screw 
dislocation in both the spatial and space-time cross- 
sections of the beam. The spatio-temporal vortex beams 
are polychromatic and can naturally appear in problems 
with non-stationary (moving) sources or media. Lo- 
calized wave-packets with spatio-temporal vortices and 
tilted intrinsic OAM can readily be constructed by con- 
sidering slightly delocalized (e.g., Gaussian) Fourier spec- 
tra instead of the delta-functions considered here. It 
should be emphasized that our analysis is valid for both 
massless and massive fields, in both relativistic and non- 
relativistic cases. 

Quite naturally, a complete relativistic family of 
spatio-temporal Bessel beams is constructed via the 
Lorentz transformations of the spatial (monochromatic) 
Bessel beams. This means that even the stationary OAM 
states of light or quantum particles will be seen as spatio- 
temporal states in the case of a transversely-moving ob- 
server or source. In this manner, the deformation of the 
phase pattern of a moving vortex [cf. Fig. 2 or 4 with 
Fig. 1] offers a sensitive interferometric tool detecting 
relativistic effects even at non-relativistic velocities. In- 
deed, purely relativistic deformations of the phase fronts 
(intimately related to the Lorentz transformation of time 
[9]) become significant for velocities v comparable with 
c 2 /ujr ~ c 2 fcj^ /u, where r ~ fcjg is the beam radius. For 
paraxial beams with 6 = k± c/uj <C 1, this allows the ob- 
servation of strong deformations at speeds v <~ 8c <C c. 
This is explained by the fact that the transverse phase 
velocity of the motion of the phase fronts is much larger 
than the speed of light, ui/k±o = 9~ 1 c 3> c. 

Spatio-temporal OAM states of light or particles can 



appear in a variety of systems involving moving frames 
and sources. In particular, they can be created upon 
scattering by moving objects, emitted by satellites doing 
quantum communications [11] or by natural astrophys- 
ical sources [13,14], and produced in collisions of high- 
energy particles [12]. Furthermore, transversely-moving 
vortex solutions can be important in two-dimensional 
wave fields, such as surface plasmon-polaritons [21]. 
While we have discussed the simplest vortex-beam config- 
urations, belonging to the Bessel-beam family, note that 
optical vortex lattices can appear from the generic in- 
terference of only three monochromatic plane waves [22] . 
In a similar manner, an infinite lattice of spatio-temporal 
vortices can appear from the generic interference of three 
plane waves with different wave vectors and frequencies. 
This offers a simple way of generating spatio-temporal 
vortices in an optical laboratory. 

In this work we only considered scalar waves and or- 
bital angular momentum. For vector (e.g., electromag- 
netic) waves carrying spin angular momentum, the rel- 
ativistic transformation properties can be more compli- 
cated. Indeed, the spin angular momentum is produced 
by a circulating spin current, which, however, does not 
transport energy [17,23]. This current is unrelated to 
phase gradients and have specific properties distinct from 
those of the regular orbital current. Although one could 
expect that the integral angular momentum of the beam 
is transformed in the same way independently of its spin 
or orbital nature, the internal deformations of the beam 
can be different. For instance, it is known that spin 
has a purely intrinsic nature, and cannot have extrin- 
sic contributions produced by the relativistic Hall effect, 
Eq. (20). Relativistic transformations of vector beams 
carrying both spin and orbital angular momenta will be 
considered in detail elsewhere. 
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